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ABSTRACT 

"<n : 

We propose that the energy source of the outburst of V838 Mon and similar 
^ ! objects is an accretion event, i.e., gravitational energy rather than thermonuclear 

runaway. We show that the merger of two main sequence stars, of masses Mi ~ 
Q ! 1.5M Q , and M 2 ~ 0.1 — O.5M can account for the luminosity, large radius, and 

! ! low effective temperture of V838 Mon and similar objects. Subsequent cooling 

CN ! and gravitational contraction lead such objects to move along the Hayashi limit, 

\ as observed. By varying the masses and types of the merging stars, and by 

considering slowly expanding, rather than hydrostatic, envelopes, this model can 

\Q ', account for a large range in luminosities and radii of such outburst events. 

■ 

o : 

Subject headings: stars: supergiants — stars: main sequence — stars: individual: 
O '■ V838 Mon, V4332 Sgr, M31 RV - stars: mass loss 

P . 1. Introduction 

£3 ■ 

V838 Mon (Munari et al. 2002, hereafter (MHK02); Kimeswenger et al. (2002)), V4332 
Sgr (Martini et al. (1999)), and most probably M31 RV (M31 Red Variable - Mould et al. 
(1990)), form a peculiar and mysterious group of erupting stellar objects. They do not fit 
to any known class of stellar outbursts. Just from the light curve they can be classified as 
slow novae but their spectral evolution clearly shows that they are not. After developing a 
K-type giant-supergiant spectrum near maximum they evolve to lower effective temperature 
and fade as very late M-type giants or supergiants. The bolometric fading at a very low 
effective temperature has undoubtly been observed for V4332 Sgr (Martini et al. (1999)) 
and M31 RV (Mould et al. (1990)). It seems that also V838 Mon is evolving in the same 
direction (e.g. Osiwala et al. (2002)). No hot phase, as e.g. nebular stage in classical novae, 
has been observed in any of the three objects. These observational facts almost certainly 
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exclude any kind of thermonuclear event models, such as runaway models for classical novae 
(see e.g. that proposed for M31 RV by Iben & Tutukov (1991)) or a late post-AGB He-shell 
flash (born- again AGB), which are usually discussed while trying to interpret the above 
variables ((MHK02); Kimesweneger et al. (2002)). The point is that the thermonuclear 
models after relatively cool phase at (visual) maximum always evolve to higher and higher 
effective temperature before they star fading in luminosity. The reason is that the stellar 
envelope, inflated by the initial thermonuclear outburst, starts shrinking well before the 
nuclear reactions start declining. 

As usually in discussions of nature of a stellar outburst the observational data on the 
pre-outburst object are very important. Unfortunantely they are very scant in our case. 
Photometric data for V838 Mon in pre-outburst are consistent with an F-type main sequence 
star ((MHK02); Kimesweneger et al. (2002)), althought observational uncertainities are 
important. For V4332 Sgr a K-type star is suggested as a progenitor (Martini et al., (1999)). 
There are no data on the pre-outburst state of M31 RV. 

If a thermonuclear event is excluded as argued above, the next energy source that comes 
to mind is an accretion (or merging) event. We therefore investigate in this paper a scenario 
in which a ~1.5 M© (F-type) main sequence star accretes a less massive objects. Given the 
energy liberated during the observed outbursts it is clear that accretion of at least a massive 
planet onto the main sequence star has to be involved. 

While constraining our model we use the results from observations of V838 Mon, as 
this is the best observed object from the three (e.g. (MHK02), Kimesweneger et al. (2002), 
Kolev et al. (2002), Osiwala et al. (2002)). This star erupted in the begining of 2002, first 
by ~6 mag. (in V) and as a K-type star it reached a luminosity of ~ 500(-D/lkpc) 2 L Q , and 
a radius of ~ 55(-D/lkpc)R Q , where D is the distance to V838 Mon. Next, after a month of 
a slow decline and cooling it had another brightning by ~4.5 mag. and reached a luminosity 
of ~ 1.0 x lO 4 (_D/lkpc) 2 L , and a radius of ~ 12O(£>/lkpc)R . Subsequent evolution 
can be described as a more or less horizontal evolution on the H-R diagram towards lower 
and lower effective temperatures somewhat interupted by a minor brightnening a month 
after the second outburst. The effective radius has countinuously been increasing and in 
April reached a value of ~ 400(-D/lkpc)R Q ((MHK02)). Numerous spectral lines have been 
observed to show P-Cyg profiles during the outburst indicating an outflow with a velocity of 
200-500 kms^ 1 , but no mass loss rate estimate is available. Lines of Lil and Ball have been 
observed, but no analysis of the element abundances has been done. In mid April the star 
started a rapid decline in V accompanied by the spectral evolution to late M-types (Osiwala 
et al. (2002)). 

A scattered light echo have been detected around V838 Mon and it was used to drive 
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its distance. While (MHK02) derive a distance of D = 0.79 kpc and Kimeswenger et al. 
((2002)) derive a distance of D = 0.66 kpc, a calcualtion by Bond et al. ((2002); H.E. Bond, 
private communication) yields a much larger distance of D ~ 3 kpc. Above, and in the rest 
of the paper, we scale the properties of V838 Mon with the distance, such that our model 
covers this entire distance scale. Indeed, one of the adavantages of our proposed model, 
is that it can account for a large range of luminosities, because the luminosity is fixed by 
the type and mass of the merging stars. Thus we can also account for M31 RV, which was 
significantly more luminous than the other two objects and reached ~ 1O 6 L 

Our paper is organized as follows. In Sect. 2 we estimate the general properties of 
a main sequence star, the primary, which accretes mass from a destructed main sequence 
companion. These results are used in Sect. 3 to proposed a scenario to the eruption of V838 
Mon. Our discussion and summary are in Sect. 4. 



2. Accreting main sequence stars 

Accreting main sequence stars were found to expand to giant dimensions in several 
different kind of calculations. These include the merger of two main sequence stars, in 
studying the formation of blue straggler stars (e.g. Sills et al. (1997)), mass accretion by a 
main sequence companion inside an envelope of a giant (Hjellming & Taam 1991, hereafter 
(HT91)), and an accreting low mass main sequence star (Prialnik & Livio 1985, hereafter 
(PL85)). 

PL85 calculate the outcome of an accreting fully convective 0.2M Q main sequence star. 
When a large fraction of the free fall energy of the accreted matter is retained, they find the 
star to expand. (PL85) parametrized the fraction of energy retain by a parameter a, 

U = aGM*M acc /R*, (1) 

where M* and /?* are the mass and redius of the accreting main sequence star, respectively, 
and M acc is the accreted mass to the envelope. For an orbiting companion that is being 
destructed and accreted, the maximum available energy has a = 0.5. Not all this energy 
will be retain by the accreted matter which build the envelope in our model. Large fraction 
of the energy goes to radiaton and mass loss observed in the eruption of V838 Mon and the 
other objects. (PL85) calculate the accretion of mass up to AM = 2.5 x 10~ 3 M Q , and find 
that when a is not too small the star expands. A good fit to their results of a — 0.1 and an 
accretion rate of M acc = 10~ 2 M Q yr _1 , and when the accreted mass is M acc > 0.5 x 1O _3 M , 
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gives for the radius of the star 



Extrapolating this fit to much higher accreted mass and radii, we find that to obtain a radius 
of ~ 50 — 15Oi? (characteristic fro the first brightening of V838 Mon), the accreted mass 
should be M acc ~ 0.02 — O.O7M . With higher accretion rate and higher values of a, (PL85) 
find larger radii. 

In the calculations of (HT91), a 1.25M main sequence star accretes from a giant up 
to M acc = 0.115M Q , at an increasing rate up to M acc ~ lO~ 2 M yr _1 . The final radius of 
the accreting star depends on the entropy of the accreted gas, being lower for lower-entropy 
accreted gas. We take the higher entropy case; we later take the low entropy-high density 
gas to be accreted into the primary star. In this case ((HT91)'s model 2c), a good fit to 
figures la and 3a of (HT91) is 

— ~ 1.4 x io 17 ( Macc / M ©). (3) 
R Q 

This is a much stronger dependence than that found by (PL85), and cannot be extrapolated 
to much larger accreted mass since then the radius becomes too large. However, this formula 
does demonstrate that a small accreted mass can cause a large expansion when the star is 
already a giant. By the last fit, to obtain a radius of R ~ 300 — 15OOi? (characteristic for 
the second, major outburst of V838 Mon), the accreted mass should be M acc ~ 0.14 — O.2M . 
The strong dependance in the last equation, suggests that with more energy supplied, as in 
the cores merger we discussed later, less mass can be inflated to large radii. Despite the 
large differences between the two works cited above, and the very different fits to (PL85) 
and (HT91), the amount of accreted mass, required to inflate a giant to R* ~ 3OOi? is 
similar, and amounts to M acc ~ O.15M . In the proposed scenario, two main sequence stars 
merge on a dynamical scale, following dynamical instability (Rasio & Shapiro (1995)) and 
the accretion rate can be as high as M acc ~ O.O3M /1 day = lOMoyr^ 1 . At such a high 
accretion rate, and with most energy deposited in the accreted mass, an accreted mass lower 
than O.15M can be inflated to giant dimensions. Therefore, for our scenario we take the 
inflated envelope mass to be M acc ~ 0.05 — O.3M . 

In the calculations cited above the accreted mass was added to the outer layers of the 
inflated star. In our scenario, the mass is added at the base of the inflated envelope, at the 
surface of the original primary star. A somewhat different treatment of the energy budget is 
required. Let the destructed companion be split into low entropy-high density gas of mass 
M c , i.e., the core of the companion, which is accreted and sinks into the core of the primary 
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star, down to a radius 77 < R\ (where R\ is the initial radius of the accreting primary star) 
and a high entropy mass M acc which goes into the inflated envelope. This assumption is 
based on the numerical calculations of Sandquist, Bolte & Hernquist ((1997)). The total 
energy of this structure is 

GM X M C 

E g ~ — + E cnv , 4 

where E cnv is the gravitational energy of the envelope. For an envelope having a density 
profile of p = Kr~ 2 - 5 from Ri to where K is a constant, which crudely fits the envelope 
found by (PL85) for the model cited above, the envelope energy has a simple approximated 
form for R* S> Ri, given by 

G(M 1+ M C )M„ /Ri\ l/2 ... 
nv ~ UJ ' <5) 

The original energy of the binary system is 

GM1M2 fR , 

E ~ -0.5 , (6) 

a 

where a is the orbital separation at destruction. In addition, an energy E 2 > should be 
supplied to destruct the companion. The condition to inflate the envelope becomes 

E Q > E g + E 2 . (7) 

If M c = 0, this condition can be met only if a > O.BR^R^/ R^ 1 / 2 . For R^/R ± ~ 500, this 
condition reads a > 10Ri. However, we do expect the destruction to occur when a < 2R\. 
We conclude that a large fraction of the destructed companion in our scenario must be 
accreted as a low entropy mass into the companion. 

As an example, consider the case of M c = kM 2) hence M acc = (1 — k)M 2 , 77 ~ 0.25i?i, 
as the size of the solar core, i?* > 300i?i, and a = 2R 1 . We also take the destruction energy 
of the companion to be 

R2 

The envelope of the companion before destruction, of mass M acc , has a smaller average radius 
than R 2 . On the other hand, there is thermal energy of the companion. Hence (3 ~ 2 — 3. 
Another energy source not considered here is the rotation energy of each of the two stars, 
which are likely to corotate in our scenario (see next section). Condition (7) becomes then 

1<164 _ 4W1 _ 4) |^ +4(1 _,)( 1 + ^)(|) 1/2 . (9) 
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From this condition it is clearly seen that when k — 0, the envelope can be inflated only to 
i?* < 1QR\. For k > 0.3 and ^> Ri, the last term, which is the energy of the inflated 
envelope can be neglected. Taking also R\M 2 ~ R2M1, the last condition is simplified to 

1 < 4fc[4 - (3(1 - k)] (10) 

for the assumed parameters here, the last condition is met by a large margin for k > 0.3. 
This means that even by relaxing some of the assumptions used here, a giant may still be 
formed from such a process. 

Not all the energy will be retain by the accreted matter which build the envelope in 
our model. Some fraction of the energy goes to the fast wind and radiation as observed 
in the eruptions. It is interesting to note that the radiation requires only small amount of 
mass to be accreted on the surface of the original primary star. Taking the energy of the 
effective-accreted mass that goes to radiation to be that of a Keplerian disk, we find this 
effective mass accretion rate to be 

For a distance to V838 Mon of D — 3 kpc, the luminosity is L ~ 1O 5 L , and the accreted 
mass over several years is still very small. 



3. The proposed scenario 

Based on the previous section, in particular on our estimate that a mass of M acc ~ 
0.05 — 0.3 is required to inflate the envelope to ~ 500 — 15OOR , we proposed the following 
general scenario for the eruption of V838 Mon. The progenitor was a close binary system 
composed of two main sequence stars: a primary with M x ~ 1.5M Q , and a companion with 
a mass 0.1M Q < M2 < O.5M . The small orbital separation implies that the spins of both 
stars were synchronized with the orbital period, and the orbit was circular. Due to evolution, 
mainly magnetic winds, possibly from both stars, the system lost angular momentum at a 
relatively fast rate, and the orbital separation decreased. One of the stars filled its Roche 
lobe, and transferred mass to the other. For main sequence stars, the increase of stellar 
radius with the stellar mass is steeper that the increase of the average Roche lobe size with 
mass (e.g., see formula by Eggleton (1983)), and the more massive primary star is likely to 
fill its Roche lobe first. At about the same time as Roche lobe overflow occurs, a binary 
system of two low mass main sequence stars becomes dynamically unstable (Rasio & Shapiro 
(1995)). This implies a very fast merger and mass transfer of the outer layers of the stars. It 
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is possible that the mass transfer of the outer layer of the primary star onto the secondary 
caused the envelope to be first inflated around the secondary star. In any case, mass transfer 
of the outer layers of the primary or of both stars, we propose, caused the first burst in 
luminosity, when V838 Mon reached a radius of ~ 5Oi? . Then, the secondary was accreted 
onto the primary, in particular the core of the secondary merged with the core of the primary, 
a process that released much more energy. We proposed that this caused the second burst 
after about a month. 

The merging process is likely to be complicated, with different epochs of variation in 
the luminosity, as parcells of gas being accreted and the merged stars, both envelope and 
the merging cores, are settled to a new equilibrium. As noted in equation (11), a very 
small amount of accreted mass can lead to a relatively strong burst. Such a variation could 
cause the third outburst, which occurred about a month after the second burst, and which 
required, if lasted one month, only ~5x 10^ 5 M Q (D/1 kpc) 2 to be accreted on the initial 
surface of the primary star. 

4. Discussion and Summary 

Our main goal is to point out that the energy source of the outburst of V838 Mon 
was, and still is, an accretion event, i.e., gravitational energy, rather than thermonuclear 
runaway. We considered the merger of two main sequence stars, of masses M 1 ~ 1.5M , 
and M 2 ~ 0.1 — O.5M . In our proposed scenario, the first outburst, with a luminosity 
of ~ 5OO(-D/lkpc) 2 L and an envelope radius of ~ 55(_D/lkpc)i? , was caused by 
the merger of the outer layers of the two stars, via dynamical instability. First mass was 
transfered from the more massive primary to the secondary; most likely the envelope was first 
inflated around the secondary. Then, after about a months, the cores of the two stars collided, 
releasing more energy, this led to the expansion of V838 Mon to i?* ~ 4OO(-D/lkpc)i? , and 
its luminosity to reache ~ 1.0 x lO 4 (-D/lkpc) 2 L . Equation (11) shows that a small 
amount of accreted mass can lead to large luminoisty variations. As the two merged stars 
settled to equilibrium, the luminoisty varies; such a process caused the third outburst as 
well. 

As mentioned in Sect. 1 a light echo have been detected around V838 Mon during 
outburst. (MHK02) and Kimeswenger et al. ((2002)) interpret this as due to scattering 
on circumstellar dust. This obviously rises a question about the origin of the circumstellar 
dust. In our scenario the dust would reside in remains of a protostellar cloud from which 
the V838 Mon binary system has been formed. This would imply a relatively young age 
of V838 Mon which is consistent with its low distance (~13 pc) from the Galactic plane 
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((MHK02)). It is, however, also possible that the echo is due to interestellar dust, as sug- 
gested by Bond et al. ((2002)), and then obviously it has nothing to do with the nature of 
V838 Mon. 

In principle, mergers of other types of stars is possible. We note that the possible 
locations of the progenitor of V838 Mon on the HR-diagram ((MHK02)), intersects the 
location of a zero-age horizontal branch (HB) star which lost most of its envelope on the 
RGB (D'Cruz et al. (1996)). Merger of such an HB star with a low mass main sequence 
star may leed to an inflated envelope around the HB star, similar in many respects to the 
scenario proposed in the previous section. The HB star lost most of its envelope because 
of the interaction with the low mass main sequence companion via a common envelope 
evolution, which started when the HB progenitor was on the upper red ginat branch. In this 
scenario the dust-halo can be attributed to the high mass loss rate on the red giant branch. 

Finally it is worth of noting that the effective photosphere of the three objects seems to 
increase more or less steadilly with time during outburst. In the case of M31 RV it expanded 
up to 8OOOR (Mould et al. (1990)). This would imply that the observed photospheric 
regions are not in hydrostatic equillibrium. The models of (HT91) and (PL85), on which 
our considerations in Sect. 2 are based, are in hydrostatic equillibrium. Slowly expanding 
envelopes would probably require less mass and thus lower M acc than hydrostatic envelopes 
inflated to large radii. Thus accretion of an object of even lower mass then estimated in 
Sect. 2 could give origin to the observed outbursts. 

This work has partly been supported by the Polish State Committee for Scientific Re- 
search through the grant No. 2.P03D.020.17. N. S. has also been supported by a grant from 
the USA-Israel Binational Science Foundation and the Israel Science Foundation. 
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